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In  1961  C.  I.  Beard  found  that  experimental  values  of  the  rough-surface  reflection 
coefficient  for  the  coherent  reflected  field  are  larger  than  the  values  given  by  the  generally 
accepted  theoretical  formula-^  A  =  exp  (-2(2irg)2]  for  values  of  g  -  (oh  sin  i‘)j\  greater 

0.1  radian.  Here  oh  is  the  standard  deviation  of  the  sea  surface  elevation,  \p  is  the  grazing 
angle,  and  X  is  the  electromagnetic  wavelength. 

The  expression  for  “K  A  is  the  Gaussian  theoretical  curve  (sometimes  called  the 
“roughness  factor”)  first  published  by  W.  S.  Ament  in  1953.  Since  then  Ament’s  result 
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20.  Abstract  Continued 

has  been  obtained  by  numerous  other  workers. 

The  disagreement  between  theory  and  experiment  has  been  unresolved  for  over  a 
decade;  and  it  has  been  agreed  generally  that  theoretical  models  based  on  geometric  optics 
could  not  yield  a  result  more  powerful  than  that  of  Ament. 

In  this  report  it  is  shown  that  the  coherent  reflected  field  is  given  by  X  *  exp 
{-2(2irg)2]/0[2(2jrg)2]  for  0  <  0.3,  where  /0(x)  is  the  modified  Bessel  function  Jq(ix). 

This  result  is  derived  using  geometric  optics,  assuming  a  spherical  wavefront  incident  on  a 
Gaussian  collection  of  sinusoidal  surface  waves  with  uniform  phase  distribution.  Further, 

X  agrees  with  Beard’s  experimental  curve,  with  a  systematic  difference  in  g  of  10%. 

Beard  has  estimated  that  wave  elevations  were  within  10%  of  their  correct  values. 

It  is  shown  how  is  obtained,  assuming  a  plane  wavefront  incident  on  a  Gaussian 
collection  of  horizontal  strips  or  equivalently  simple  functons. 
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GEOMETRIC-OPTICS  THEORY  FOR  COHERENT  SCATTERING  OF  MICROWAVES 

FROM  THE  OCEAN  SURFACE 


INTRODUCTION 

In_1961  C.  I.  Beard  (1)  found  that  experimental  values  of  the  coherent  reflected 
field  \ElEfr  T|  are  larger  than  the  values  given  by  the  generally  accepted  theoretical  formula 

E  \  9 

J^\$A  =  exp  [-2(27rg)2J 

for  values  of  g  =  (oh  sin  i//)/X  greater  than  0.1  radian.  Here  £  is  the  average  electric 
field  due  to  the  “sea  surface”  $A,  E&  is  the  field  due  to  the  direct  wave,  V  is  the  smooth 
sea  reflection  coefficient,  oh  is  the  standard  deviation  of  the  sea-surface  elevation,  i//  is  the 
grazing  angle,  and  X  is  the  electromagnetic  wavelength  (Fig.  1). 


R 


The  expression  is  the  Gaussian  theoretical  curve  (sometimes  called  the 

“roughness  factor”)  first  published  by  W.  S.  Ament  (2)  in  1953.  (Ament  claims  that 
the  result  was  derived  by  Pekeris  and,  independently,  by  MacFarlane  in  the  40’s.)  Since 
then  Ament’s  result  has  been  obtained  by  numerous  other  workers.  (See  for  example 
Ref.  3.) 


Manuscript  submitted  December  14,  1973. 
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The  disagreement  between  theory  and  experiment  has  been  unresolved  for  over  a 
decade;  and  it  has  been  agreed  generally  that  theoretical  models  based  on  geometric 
optics  could  not  yield  a  result  more  powerful  than  that  of  Ament. 

In  this  report  we  show  that  the  coherent  reflected  field  is  given  by 

~  =  exp  l-2(2ng)2]I0[2(2wg)z] 

E6r  $ 

for  0  g  <  0.3,  where  /0(x)  is  the  modified  Bessel  function  J0(bc).  We  derive  this 
result  using  geometric  optics,  assuming  a  spherical  wavefront  incident  on  a  Gaussian 
collection  S  of  sinusoidal  surface  waves.  Further,  |£/£§r|S  agrees  with  Beard’s  experi¬ 
mental  curve  (Fig.  2),  with  a  systematic  difference  in  g  of  10%.  Beard  has  estimated  that 
wave  elevations  were  within  10%  of  their  correct  values. 

We  show  how  \EIEsV\Sa  is  obtained,  assuming  a  plane  wave  front  incident  on  a 
Caussian  collection  of  horizontal  strips. 


BASIC  ASSUMPTIONS  AND  DEFINITIONS 

Suppose  that  T  is  an  isotropic  emitter  of  electromagnetic  waves  with  wavelength  X 
situated  at  a  distance  h2  above  file  mean  sea  surface  and  that  A  is  an  isotropic  receiver 
which  is  at  a  distance  /ij  above  the  mean  sea  surface.  We  take  the  mean  sea  height  to  be 
equal  to  zero.  Suppose  that  the  horizontal  distance  between  T  and  R  is  d,  where  d  is 
small  enough  that  the  earth’s  curvature  need  not  be  considered  (Fig.  3).  We  assume 


Fig.  2— Comparison  of  theoretical  and  experimental  results 
for  coherent  forward  scattering 


that  the  sea  surface  can  be  described  by  some  collection  of  functions  CS  .  To  be  more 
specific,  we  assume  that  at  any  given  moment  the  functions  which  we  need  to  describe 
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the  sea  surface  locally  are  in  CS  .  Let  T'  be  the  image  of  T  with  respect  to  the  tangent 
line  through  P,  a  typical  reflection  point  on  the  sea  surface.  Let  Xr?2  -  kj)  be  the 
coordinates  of  T',  so  that  the  coordinates  of  Pf  are  (X^/2,  Xi?2/2).  Note  that  when  P  is 
on  the  s  axis  and  we  have  specular  reflection,  rjx  =  rj2  =  0,  i.e.,  the  image  of  T  through 
the  s  axis  is  at  T's(0,  -h2).  When  we  have  this  situation,  let  ^  be  the  angle  TPO.  It  is 
easy  to  see  that  d  tan  $  =  hx  +  h2.  We  remark  that  we  will  be  using  two  independent 
coordinate  systems,  viz.,  the  s,  y  system  and  the  rjlf  rj2  system  with  origin  T's.  Note  that 
the  coordinates  of  O  in  the  s,  y  system  are  (0,  h2)  in  the  i?j ,  rj2  system. 

In  this  work  we  use  geometric  optics,  that  is,  we  assume  that  only  waves  reflected 
from  a  favorable  slope  will  reach  the  receiver  and  that  the  local  angle  of  incidence  will 
equal  the  local  angle  of  reflection.  Further,  we  assume  that  all  reflections  occur  in  the 
s,  y  plane.  We  also  will  assume  that  the  sea  surface  is  “slightly  rough,”  so  that  over  a 
long  period  of  time  the  collection  of  points  T',  the  image  of  T  =  il  -p ,  will  form  a  closed 
connected  set  which  has  no  holes,  that  is,  a  closed,  simply  connected  region.  It  will  also 
be  assumed  that  T'  is  an  interior  point  of  .1  -p . 

The  following  definitions  will  be  helpful: 

=  /tj  +  ^2  '  ^2 
*2  =  2h2  -  Xr*2 
*3  =  d  ~  H 
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if  -  d2  +  (/i,  +  A2)2 


yl  -  -2ir  cos 
y2  =  -2ir  sin  \J/ 

Note  that  rt  is  the  specular  distance  RT*t.  The  distance  RI*  in  general  is  r  -  (£2  +  £§)1/2. 

We  will  show  later  that  under  suitable  conditions  7^  is  the  center  of  symmetry  of  <jj 
and  that  if  (Xrjlt  Xr?2)  T,  then 


Xr?j 

<  1, 

Xt?2 

d 

h2 

so  that  J  T  is  bounded  away  from  T  and  R.  Let  be  the  smallest  rectangle  which  con¬ 
tains  i  7  and  whose  sides  are  parallel  to  the  coordinate  axes.  Then  in  view  of  the  previous 
inequalities,^  is  a  neighborhood  of  T't,  i.e.,  the  dimensions  of  fR  are  small  compared  to 
d  and  h2.  We  then  are  justified  in  expanding  r  in  a  Taylor  series  about  tj1  =  rj2  =  0  and 
dropping  terms  higher  than  those  of  first  order  in  rjx  and  tj2,  giving 


2irr 

—  “  yo  +  W  +  W  • 


(1) 


If  we  define  the  row  vector  y  =  [7lf  y2] 


and  let  rj  = 


,  we  can  write  Eq.  (1)  as 


2ir  r 

T  *  7°  +  y11  • 


(2) 


CALCULATION  OF  THE  ELECTRIC  FIELD  DUE  TO  A  SINGLE  REFLECTED  RAY 

The  electric  field  E  due  to  a  single  reflected  ray  is  given  by 

2irr 
—  + 

X 

where  T  is  the  complex  reflection  coefficient  and  £§  is  the  electric  field  due  to  the  direct 
wave  which  travels  path  length  5  given  by 


/  27r5 


E  =  £5  —  |  f  |  exp  ij  exp 


arg  T  i 


(3) 


62  =  d2  +  (h2  -  hx)2  . 


4 


w.*as.^pwjWU|ir»g4^M.jj * ■-^  <aM»y^i^.^»wwgfgwPCT^^TO1^gTwW3^^^M^^wp^".f?»w.JTO^BTyg-i ; g » 


NRL  REPORT  7705 


We  consider  now  the  quantity  6/r  and  state  the  following. 
Observation:  If  2hj/d  <  1  (<  means  leu  than  or  of  the  order  of). 


and 


then 


h\ 

T  *  *• 


s 

—  ~  1 
r 


Proof:  Using  Eq.  (1), 

1  =  _ [d2  +  (h2  -  hy  )2  j  1/2 

r  fd2  +  (h2  +  )2ja/2  _  Xt?j  cog  ^  5^  ^ 


\*  +  (*2  -  M2' 

1/2 

\ 

?-+  +  M2 

1  - 


Xt?, 


[d*  +  (A,  +  Aj)*]i/*  C0§^  •  jg"T(^  +\)2]i/2  sin* 


Now 


Xrjj 


Xt?j 


and 


l(«r*  +  ( h2  +  /tjjzjmj  v  |- 

d  1 

1 

X»?2 

Xl?2 

Xr?2 

(d2  +  (h2  +  hj)2]i/2 

s 

h2  +  hj 

< 

s 

«  l; 


therefore 


—  %  d2  +  +  *2  -  2^1^2 
^  *“  d2  +  h2  +  *2  +  2hjh2  • 

However,  2h2/d  <  1  and  hx  Id  <  1  imply  that  2hlh2/d2  <1,  so  that  5/r  «  1. 
We  then  are  justified  in  writing  Eq.  (3)  approximately  as 


E  ~  Es  |rl  exp  (  ~-ij  exp 


2  nr 


*  **x  r )  i 


(4) 


5 


BROWN  AND  MILLER 


It  is  not  hard  to  see  that  T  =  HtJj,  rj2).  Since  $  x  C  and  !R  is  small,  we  will  make  the 
simplifying  assumption  that  T  =  T(0t  0)  =  r(\p)  in  i  j.  For  the  surfaces  we  consider  later 
the  point  (0,  0)  will  be  the  center  of  symmetry  of  i  7,  so  that  T  =  T(^)  can  be  viewed  as 
some  average  in  ^-p.  Combining  Eqs.  (2)  and  (3)  and  defining 


£5ir(i//)|exp 


exp(-70i)exp  [-iarg  T(^)j 


gives 


E  =  Eq  exp  (-iyn)  •  (5) 

We  remark  that  2h2/d  <  1  and  h^/d  <  1  imply  that  (h^  +  h2)/d  <  1/2,  i.e.,  tan  ^  <  1/2, 
so  that  <  21°.  Hence  Eq.  (5)  is  valid  provided  that 


and 


1, 


<  1, 


\p  <21°  . 


A  RESULT  DUE  TO  AMENT 

Suppose  $  A  is  the  class  of  functions  y  =  H,  where  the  sea  elevation  H  is  distributed 
normally  with  standard  deviation  aH  and  \2H\  <  min  {hly  h2}.  This  is  equivalent  to  the 
class  of  functions  Ament  used  in  Ref.  2,  which  is  the  class  H(s)  where  H  “varies  so  slowly 
with  s  that,  in  the  neighborhood  of  each  s,  the  surface  is  approximately  a  plane  parallel 
to  the  s  axis.”*  It  is  easy  to  see  that  for  the  class  S  A, 


•An  equivalent  way  of  defining  the  class  is  the  following: 

Let  U  }  be  a  disjoint  collection  of  open  intervals  such  that 

U',  -  «-S. 

/»i 

where  S  is  a  set  of  measure  zero  and  A  is  the  real  line.  Define 

r  1  if  X  e  / 

K<X>  = 

(.0  if  X  «  / 

Let  j  be  a  sequence  of  real  numbers,  and  define  a  simple  function 

AX)  ■  L]  «K(X). 

/■1 


Let be  the  class  of  all  simple  functions  (see  Fig.  4). 
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*T  *  j(Xl»1,Xl?2)!t?i=0tl?2!5  ~  f  . 

so  that  in  this  case  =  d'p  a  degenerate  rectangle.  Hence  7 n  “  72^2  =  -4*#  (sin  ^)/X, 
so  that  from  Eq.  (5)  we  have 


/ 4niH  sin 

E(H)  =  Eq  exp  ^ - - - j 

„  (4irH  sin  \l/\  .  (4itH  ain  4>\ 

*  — r 


The  average  electric  field  E  due  to  the  class  is  given  by 


If  we  substitute  Eq.  (6)  into  this  equation,  it  is  easy  to  see  that  we  get 


To  evaluate  these  integrals,  we  will  let  H0  -*•  +«> .  The  error  we  make  will  not  be  large, 
since  the  cosine  function  is  bounded  and  the  distribution  tends  to  Kero  quickly  when  aH 
is  small.  Performing  the  required  integrations  and  noting  that  |E0 1  =  I  £5  P(i//)|  gives 


Fig.  4  —Typical  simple  function 
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E 


h? 


(7) 


Ament  derived  this  result  for  a  plane  incident  wave  and  the  class  CS  A .  We  have  made  the 
assumption  that  the  incident  wave  front  is  spherical.  However,  if  we  are  considering  a 
plane  wave,  then  f  ==  T(^)  in  ,  and  we  have  Ament’s  result. 


CALCULATION  OF  THE  IMAGE  DUE  TO  THE  CLASS  $ 

We  assume  now  that  the  sea  surface  can  be  described  by  a  Gaussian  collection  S  of 
sinusoidal  waves  y  =  H  sin  [(2it/n)(s  +  s  )] ,  where  the  sea-wave  peak  height  H  >  0  is 
normally  distributed,  p  is  the  sea-wave  wavelength,  and  s'  is  the  sea-wave  phase,  which  is 
distributed  uniformly. 

Let  (3=  (27r//u)(s  +  s  ),  and  note  (Fig.  5)  that  the  coordinates  of  P  are  (s,  H  sin  0). 

The  image  j  j  is  determined  by  the  following  three  equations: 

Sj  (d  -  s)  =  £3(fy  -  H  sin  0);  (8) 

27rH£2  cos  0  =  ;  (9) 

2rrf/(XT/,  -  2s)  cos  0  -  M(Xi72  -  2H  sin  0) .  (10) 


Equation  (8)  is  the  condition  that  T',  P,  and  R  lie  in  a  straight  line;  Eq.  (9)  is  the  con¬ 
dition  that  the  line  through  P'  and  P  is  perpendicular  to  the  line  through  T  and  T\  and 
Eq  1 10)  is  the  condition  that  the  slope  of  y  =  H  sin  0  at  P  is  equal  to  the  slope  of  the 
linf  through  P'  and  P. 

We  now  state  the  following  observation. 
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Observation:  If 
4ir  sH 


„  2 H  2h0  h,  +  h0 

<  1 ,  —  <lt  ~1<1,  -I— £  <  l , 


d  ~ 


(11) 


then  i  -p  is  approximately  an  ellipse  and 


_  f  2H  2nh2  2H  [  2nh0  ,  \  ] 

oi»T  =  |(Xi7i,Xt?2)It?1=  cos/3,i?2  =  —  |sm/3 - cos  \p  cos /3j  |  , 

where  /?  G  [0,  2?r] ,  (3  means  "the  boundary  of. ") 

We  remark  that  the  relations  (11)  aie  equivalent  to  the  following: 


nH  2h2  h.  +  h2  ^  ^ 

2tt—  <1,  <  1,  <  1, 

/i  d  d 


(12a) 


K  = 


2  ( 2ff —  cos  (3 


2  4tts//  2// 

cos  p  +  - —  sin  0 

2 


h2p  '  hn 


<  1,  (3  G  [0,  2n]  .  (12b) 


The  relations  of  (11)  imply  those  of  (12)  since  3s  3  s/h2  >  1  (for  example,  when 
H  =  0,  s/h2  =  d/(/i1  +  h2)  >  1,  so  that  certainly  y/2s/h2  >~1,  and  2 s/h2  >  1;  thus 


y/2  1 »  2?r-^  <  1 .  K  <  1 


Relations  (12)  imply  (11),  for  when  (3  =  7r/2,  we  have  2 H/h2  <  1. 
When  (3  =  0, 


Hence  we  must  have 


AirsH 

h2p 


<  1 . 


or 


2K)2  «*• 


4ttsH 


J  H\2 

AnsH 

2  2tt— 

\  v  1 

h2p 

which  implies  that 

H 

s 

2jt —  ~ 

— 

P 

^2 
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2  it—  <  1 , 


and  so 


T-  <1. 


a  result  which  is  not  always  true.  Hence 

4  n$H 
h2n 


<  1. 


Further,  if  we  eliminate  the  parameter  /3,  we  have 


33t  -  |  (Xi?i ,  X??2)| 


where 


A  = 


+  cot2  \p  cot  \p 


2' 

L  cot  ^  1 

Proof:  From  Eqs.  (9)  and  (10)  respectively, 

„  H  B  Xtl  1 

2jt—  COS  (5  =  — ^  - - - r  , 

^  h2  l2_ 


and 


Xt?2  2H 


H 


-  sin/J 


2rr-  cos/3  =  J3L 

M  Xrjj  2s 

^2  ^2 

Eliminating  Xrjj  /h2  from  these  equations  gives 

Xt?2  2H  . 


J] 

2rr —  cos/3  = 


*2  -^SM 


v(s 


Xr?5 

hn 


2?r —  [  2 - —  )  cos  /3  — — 


2s 


(13) 
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(14) 


(15) 


(16) 


(17) 
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- (- 7VV 


d .  (18) 


+  h2 


Clearly 


Xr?, 


*1  +  h 


<  1, 


and  since  from  Eq.  (14) 


Xr?! 


„  //  2A2 

2rr - -  cos  /J 

H  d 


<  1, 


we  can  write  Eq.  (18)  approximately  as 

d 


8  = 


~ sin/3  -  Aj)  +  d 

Ai  +  A2 

*  cot  ^  (A2  +  H  sin  0) . 

From  Eqs.  (17)  and  (19),  we  have 

H 


-  H  sin/J  +  2rr 


U  \2  I  H 
—  cos  0J  ft  2  -  f  27T— -  co8/3  )  h2  cot  \J/ 


-  ^2rr-^-  cos  pj  H  cot  ^  sin  0 . 


Since 


2rr 


H  \ 

—  COS0j  // 


cot  ^  sin  0 


/  // 

[2rr—  cos|3  }  A„  cot  ^ 

\  M  1  l 


H  sin  0 


<<  1 


(19) 


(20) 


and 
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we  can  write  Eq.  (20)  approximately  as 


Xr?2 

2 


H  sin  0  - 


h2  cot  , 


and  so 


** 2 


cot  \p  COS  0 


')• 


We  easily  can  verify  that 


rj'Ari  = 


o 


(21) 


and  check  that  it  is  indeed  an  ellipse.  The  family  of  ellipses  is  shown  in  Fig.  6.  If  we 
let0  -  jt/2,  3rr/2,  we  obtain  the  interesting  ^result  that  il-p  (Class  &  )  O  {y  axis}  =  4^ 

(Class  ),  so  thct  the  major  axis  of  the  ellipse  (2///A)2  =  tj'Ajj  must  be  greater  than  4///A. 


AVERAGE  ELECTRIC  FIELD  DUE  TO  THE  CLASS  S 

For  a  given  H  and  0,  the  electric  field  due  to  a  single  reflected  ray  is  given  by  Eq. 
(5):  E  =  E0  exp(-«7Tj).  We  easily  find  that 


4ir//  .  ,  .  a 
Ttj  =  — - —  sm  w  sin  0  , 

A 


so  that 


E(H,  0)  =  E.  exp  { -  tOJL  sin  i p  sin  0 


)• 


(22) 


When  0  =  rr/2,  3?r/2,  we  get  Eq.  (6).  The  physical  interpretation  of  yi}  is  shown  in 
Fig.  7. 
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We  now  define 


=  — 


r0~2n 

E(H,  0)  di 

4o _ 

p/3«2ir 

I  di 


Equation  (23)  gives  the  average  elective  field  due  to  all  reflections  such  that  T'  is  in 


We  recall  that 


so  that 


2ir 

p  =  —  (« ♦  n 


s  =  cot  ^  (h2  +  f/sin  0) , 


—  0  *  /j0  cot  4/  +  H  cot  \p  sin  0  +  ? 
2jt  z 


Taking  differentials  in  this  last  equation  gives 


— d0  =  H  cot  ^  cos/3  dp  +  ds  , 

2ff 


so  that 


\2ff 


cot  i//  cos  /3 


Hence 


£(H)  = 


cot  ^  cos/3  )d/3 


I  -  //  cot  cos  p\ 

4o  V2tt  / 


1 


f*3’2"  „  „  //cot^  r^*2* 

£(H,  /3)  d/3 - —  El 

'0*0  ^  ^3-0 


rP“2n 

E(H,  P)  cos  p  dp  . 

Jft-n 
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But  E(H,  p)  is  a  function  of  sin  P,  so  that 

pP*2n  rP"2n 

E(H,p)co$pdp  =  E(H,  p)  d  sin  P  =  0 . 

*7?»o  iJ«o 

Hence 

1  r0*2?r 

S(H)=—  E(H,p)dP , 

Zff  J0.o 

and  on  performing  the  required  integration,  we  obtain 

£(//)  =  Vo (^“  sin*)  *  (24) 

where  Jq  >s  the  Bessel  function  of  order  0.  Note  that  n  does  not  appear  in  Eq.  (24). 

Now  let  E  be  the  average  electric  field  due  to  the  Gaussian  collection  of  all  sinusoidal 
weves,  i.e., 


where  2 H0  <  min  {hx,  h2  }  and  aH  is  the  standard  deviation  in  the  peak  height  H. 
Substituting  Eq.  (24)  into  Eq.  (25)  gives 


As  before,  we  let  H0  ~+  +°°  without  making  a  large  error,  and  on  performing  the  required 
integrations  (and  noting  that  \E0  \  =  lE^Tt),  we  obtain 


E 

s  pv  n 

/  2i wH  sin  p  \2 

/ 

( 2 naH  sin 

CAp 

LI  x  /  J 

[\  x  /  J 

where  /0(*)  is  the  modified  Bessel  function  <l0(ix). 
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If  ah  is  the  standard  deviation  of  the  sea  elevation,  then,  since  we  have  considered 
sinusoidal  waves,  it  is  clear  that  aH  *  (Fig.  8).  Hence  we  can  write 


E 


=  exp 


j2*oksin+y 

r 

J2*ch  sin  *\2 

r2l  x  l 

70 

r\  x  /J 

Fig.  8— Relationihip  between  standard  deviation  in 
peak  height,  aH ,  and  standard  deviation  in  elevation, 
oh,  for  a  sinusoidal  wave 


We  remark  that  the  use  of  sinusoidal  surface  waves  is  a  mechanism  for  obtaining  a  rich 
source  of  wave-height  wave-slope  pairs.  Further,  it  seems  that  we  obtain  enough  of  these 
pairs  distributed  in  the  right  way  to  describe  the  sea  surface  over  a  long  period  of  time. 

To  sum  up  we  provide  the  following  tabulation: 


Table  1 

Coherent  Field  for  Various  Sea-Surface/Incident-Wave  Combinations 


Sea  Surface 

Incident  Wave 

Result 

Simple  functions 

Plane 

Ament:  \E/E^r\^A 

Simple  functions 

Spherical 

Ament:  \EIE^r\^A 

Sinusoidal  waves 

Plane 

Ament-like:  \E/E^r\^ 

Sinusoidal  waves 

Spherical 

\E/E^r\h 

The  Ament-like  result  is  obtained  by  substituting  oH  =  \f$oh  in 


giving 


E 

■ 

/2ir a„  sin  \p\2 

=  exp 

-2 

\  X  )  J 

E 

=  exp 

-4 

/2 irah  sin  i^\2 

cS; 

- 

\  X  /  . 
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This  is  justified  since  when  the  incident  wave  is  plane,  the  only  reflections  that  can  occur 
are  off  the  peaks  and  troughs  (assuming  no  shadowing)  of  the  waves  (Fig.  9). 


A  »\X» 

/\y\y\7 


Fig.  9— Reflection  point*  for  a  plan  incident  wave 


A  RESULT  DUE  TO  BECKMANN 


Since  the  completion  of  this  work,  it  has  been  realized  that  another,  but  different, 
theoretical  treatment  based  on  geometric  optics  had  been  done  previously  by  Beckmann 
(4)*.  Beckmann’s  report  was  undertaken  to  try  to  calculate  theoretically  the  effect  of 
spherical-wavefront  illumination  for  comparison  with  Beard’s  experimental  results  of  the 
statistics  of  the  incoherent  field  scattered  from  a  random  water  surface  (5).  In  the 
process  of  calculating  the  incoherent  Reid,  Beckmann  also  derived  expressions  for  the 
coherent  field.  It  seems  that  these  were  not  pursued  at  the  time  since  the  interest  was 
in  the  statistics  of  the  incoherent  field.  Beckmann  obtained 


*  K  „/l  3  \ 


where  lFl(  1/2;  3/2;  K)  is  a  confluent  hypergeometric  function.  We  have  extracted  this 
result  from  Ref.  4,  Eqs.  (49),  (50),  (51),  and  (85).  Numerical  computations  indicate  that 


E 

< 

E 

< 

E 

es  r 

A 

Et,  r 

0 

for  0.1  <  g  «  0.3  . 
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